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ABSTRACT 

Proposed scaling relations of a characteristic timescale in the X-ray power 
spectral density of galactic and supermassive black holes have been used to ar- 
gue that the accretion process is the same for small and large black holes. Here, 
we report on the discovery of this timescale in the near-infrared radiation of 
Sgr A*, the 4 ■ 10^ Mq black hole at the center of our Galaxy, which is the most 
extreme sub-Eddington source accessible to observations. Previous simultaneous 
monitoring campaigns established a correspondence between the X-ray and near- 
infrared regime and thus the variability timescales are likely identical for the two 
wavelengths. We combined Keck and VLT data sets to achieve the necessary 
dense temporal coverage, and a time baseline of four years allows for a broad 
temporal frequency range. Comparison with Monte Carlo simulations is used to 
account for the irregular sampling. We find a timescale at 154lgy^min (errors 
mark the 90% confidence limits) which is inconsistent with a recently proposed 
scaling relation that uses bolometric luminosity and black hole mass as parame- 
ters. However, our result fits the expected value if only linear scaling with black 
hole mass is assumed. We suggest that the luminosity-mass-timescale relation 
applies only to black hole systems in the soft state. In the hard state, which is 
characterized by lower luminosities and accretion rates, there is just linear mass 
scahng, linking Sgr A* to hard state stellar mass black holes. 

Subject headings: black hole physics. Galaxy: center 
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Introduction 



Cosmic black tioles (BH) sliow a very wide range of masses: from stellar masses to 
hundreds of millions of solar masses. An open question has long been whether the accretion 
and variability processes occurring in the immediate vicinity of the event horizon are the 
same over this broad mass range. To test this intriguing possibility, investigators identified 
and studied a characteristic timescale associated with the aperiodic X-ray variability of black 



hole X-ray binaries (BHXRBs) and active galactic nuclei (AGN; see, e.g., Uttley et al. 2002 
Markowitz et al. 2003 [Uttley fc McHardy 2005). This timescale corresponds to a break in 
the power spectral density (PSD) at a certain temporal frequency where a power-law of slope 
7 (with P{f) oc /"''') breaks to a steeper slope /3 > 7. 



McHardy et al.| (2006) proposed a scaling relationship between the break frequency, the 
mass of the black hole and the bolometric luminosity of its accretion flow. This extended 
earlier work which hypothesized that the break timescales of AGN scale linearly with BH 
mass from the timescales observed in BHXRBs, albeit with some scatter. [McHardy et al. 



(2006) found that this scatter can be explained by introducing the bolometric luminosity 
as a correction factor. This lead them to conclude that AGN are scaled-up galactic BHs. 
However, their sample consists of only 10 AGN and it is therefore desirable to test their 
scaling relation with newly determined AGN break frequencies. The BH in our own Galactic 
center (identified with the radio source Sgr A*) is an an especially interesting object, as 
it is the most underluminous BH accretion system observed thus far (with a bolometric 
luminosity nine orders of magnitude lower than its Eddington luminosity) . It can therefore 



test the McHardy et al. (2006) relation in so far inaccessible regions of the parameter space. 



The existence of a supermassive BH in the center of the Milky Way has been demon- 
strated beyond reasonable doubt by the proper motion of stars detected in the near-infrared 



(NIK) waveband ( 


Eckart & Genzel 


1996 


Ghez et al. 


1998 


2000 


2000 [Schodel et al. 


2002 


Gillessen et al. 


2008 


). X-ray and NIR 



Sgr A* has been observable since 2001 and 2003, respectively (Baganoff et al. 2001 Genzel 



et al. 2003 Ghez et al. 2004), which showed that it is highly variable at both wavelengths 



Eckart et al. 


2006a|b 


Meyer et al. 


2006a|b 


Trippe et al. 


2007 



X-ray monitoring campaigns revealed that each X-ray flare is accompanied by a NIR flare 
with zero time lag, which leads to the conclusion that the X-ray photons are being produced 



al. 2004 2006a 2008; Marrone et al. 2008 Yusef-Zadeh et al. 2006 Hornstein et al. 2007). 



in the NIR ( 


Eckart et 


Hornstein et al. 


2007). 



The fact that the same bunch of electrons is the source for both the NIR and X-ray flux 
makes it possible to interpret the NIR timing studies reported here in the context of the 
AGN and BHXRB X-ray results. In fact, since the X-ray background is very high due to 
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the lack of resolving power of X-ray telescopes (and thus only during the brightest flares is 
Sgr A* actually observed; see, e.g., Baganoff et al.| 2001 Belanger et al. 2005), a detailed 
study of the PSD is more complicated in the X-rays than in the NIFQ 

Recent studies of the NIR properties of Sgr A* showed that the variability on timescales 
of minutes to hours is completely described by an (unbroken) power-law PSD with a slope 
of -1.6 to -2.5 ( [Meyer et aL][2008| |Do et al.|[2508[ see also Fig. 18 in Eckart et al. 2006a.). 
In this Letter, we extend the time baseline up to years by combining Keck and VLT data. 
We report for the first time the existence of a power-law break frequency in the NIR PSD of 
Sgr A*, and we show that its deduced range is inconsistent with the McHardy et al. (2006) 
scahng relation, but fits the expected value when linear scaling with mass is assumed. 



2. The data 

The NIR adaptive optics instruments at Keck II (NIRC2) and at the VLT UT4 (NACO) 
have been used for Sgr A* observations since 2003. A time baseline of years and a sampling 
timescale of minutes are needed to analyze the broadband PSD of Sgr A* from minutes to 
years. Obviously, large gaps in the sampling pattern are unavoidable as the NIR observations 
have to be carried out during the night and available telescope time limits realistic data sets 
to a couple of nights each year. 

To cover the high-frequency to mid- frequency part of the PSD (meaning timescales of 
minutes to days in this context), we looked for Keck and VLT data from consecutive nights 
with individual night observations that are as long as possible. We identified the nights 
from 2004 July 06 - 08 (all times UT) as the most suitable. Beginning on July 06, Sgr A* 
was observed from 07:50 to 10:35 UT with Keck II, from 23:19 to 04:16 (July 07) with the 
VLT, then again from 06:35 to 10:30 with Keck II, and finally from 00:53 (July 08) to 06:53 
with the VLT (the Keck and VLT data have been published in Hornsteiii||2007 Eckart et al. 



2006a, respectively). Clearly, this dense coverage is only possible by combining the Keck and 
VLT data sets. Please note that these alternating observations were coincidental and not 
arranged. As the Keck data for these two nights were taken at L' (centered at 3.8 /xm) and all 



other data were taken at K-band (centered on 2. 15 /im)[^ we adopt the finding of Hornstein 



^We want to note that also in the NIR there might be a background due to unresolved stellar sources at 



the position of Sgr A*, which has been estimated to contribute up to 30% to the observed flux when Sgr A" 



is at its lowest flux levels (Do et al. 20081. This low background level is negligible for the purpose of our 



study. 



^More precisely, Keck uses a K' fllter (A = 2.12 /xm) and VLT a Kg fllter (A = 2.18 /zm). We caUbrated 



et al. (2007) that the NIR flux is consistent with a constant spectral index of a = —0.6 
{Fu oc z/"). Thus, we scaled the flux at L' down to the K-band level with this relation. The 
resulting lightcurve for the 3 nights is shown in Fig. [Tj All fluxes are de-reddened using 
Ak = 3.2 as it was done in Hornstein et al. (2007). Other authors have claimed variable 



and/or somewhat steeper spectral indices (e.g. Gillessen et al. 2006). However, our results 
presented here are not very sensitive to the exact value of a. 

For the coverage of the mid-frequency to low-frequency part of the temporal spectrum 
(timescales of days to years) we averaged the data of individual nights. This leads to a 
lightcurve with irregular sampling of down to one day that covers a time baseline of 4 years. 
The details of the whole data set are given in Table 1. 

The data reduction was standard, i.e. flat fielding, sky subtraction, and correction for 
bad/hot pixels. The VLT data have been deconvolved with point-spread functions (PSF) 



extracted from the individual images (Diolaiti et al. 2000). Aperture photometry was done 



on each image and the flux was calibrated relative to sources in the field with known flux. See 



Meyer et al. (2008) for more details. For the Keck data the individual PSFs have been used 



to fit sources in the field; see Do et al. (2008) for details. The two photometric techniques 



lead to indistinguishable results as shown in Meyer et al. (2008). 



The irregular sampling with the large gaps makes standard Fourier transform techniques 
unsuitable for our data set. We therefore employ the first order structure function (SF), 
defined as V{t) =< [s(t + r) — s(t)]^ > with s(t) being a measurement at time t, to look for 
a break in the power-law PSD (a power- law in the PSD translates into a power-law in the 
SF albeit with a different slope, see, e.g., [Simonetti et al.|[T985[ [Hughes et al^|1992[[Do et al. 



2008 ). The SF for our complete data set is shown in Fig. 2 The diamonds represent the high 



to mid-frequency data shown in Fig. [T} and the crosses show the mid- to low- frequency part. 
The overlapping region was used to normalize the low-frequency part such that a continuous 
SF emerge^ It shows the superposition of a constant at short lags (< 1 min) which is the 
result of uncorrelated measurement noise, a power-law portion between lags of minutes to 
hours, and a plateau/shallower power-law for long time lags (> 1000 min). This plateau at 
low frequencies points to a break in the power-law at shorter lags. To quantify the evidence 



both data sets in the same way and thus the difference is negligible. 

■^Equivalently, we could renormalize the high-frequency part. The reason why the SF is not automatically 
continuous goes back to the fact that we average the data of a whole night for the low-frequency part (which 
reduces the variance), whereas the high-frequency data are not averaged. Contrary to the periodogram 
(Markowitz et al. 20031, the SF does not have any normalization factor which accounts for that, and thus 
we use the overlapping region to determine the renormalization factor. 
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for the existence of a break frequency and its value, Monte Carlo simulations are required, 
as the sampling pattern severely distorts any intrinsic power-law PSD, and introduces the 
effects of red-noise leak and aliasing (e.g., Uttley et al. 2002). 



The simulations 



We used the following approach for our MC simulations, which is very similar to the 
PSRESP method byfUttley et al.|(|2002D and|Markowitz et al.|([2003|: 1. An intrinsic (broken 



or unbroken) power-law PSD model is assumed and corresponding lightcurves are generated 



using the algorithm by Timmer & Konig (1995). This is done independently for the high 



frequency and low- frequency part. These lightcurves, which are significantly longer and more 
densely sampled to account for aliasing and leakage, are subsequently resampled according 
to the sampling function of the real data and uncorrelated measurement noise is added (see 



Do et al. 2008, for more details). 2. The SF from the simulated lightcurves are calculated in 
the same way as is done for the data (the renormalization is done for each pair of long and 
short-term simulations individually). 3. An observed value is determined by using the 
observed SF, the mean simulated SF, and error bars equal to the rms spread of the individual 



simulated SFs at each time lag (Uttley et al. 2002). 4. The goodness of fit is computed by 



modehng the distribution for each assumed PSD model, i.e. several thousand values 
are calculated using the individual realizations of the simulations instead of the observed 
data. The probability that the model PSD can be rejected is then given by the percentile 
of the simulated values exceeded by the value of observed x^- 5. The steps above are 
repeated to scan a range of break frequencies and power-law slopes. 

For each model PSD we used 100 simulations for each the high- and low-frequency part. 
By combining both sets we arrived at 10,000 simulations to determine the x^-distribution. 
The class of models we employed as the intrinsic PSD is the singly broken power-law: 



A{f/hr)-\ f<hr, 
AU/hr)-^ f>hr. 



Here, fbr is the break frequency (the characteristic timescale), A the PSD amplitude at fbr, 
j3 the high-frequency power-law slope, and 7 is the low-frequency power-law slope with the 
constraint 7 < /5. We tested a (3 range of 0.5 — 3 in increments of 0.1, a 7 range of — 1 
in increments of 0.1, and a break frequency range of 10~^ — 3 ■ lO^^min"^ in multiplicative 
steps of 1.2. We also tested unbroken power-law models {P{f) oc /~^) which are, however, 
rejected with a likelihood of acceptance of 0% (i.e. none of the 10,000 simulations equalled 
or exceeded the observed x^)- 
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Results &6 Discussion 



The best fit result with a likelihood of acceptance of 



-+5 
-2.9 



70 can be found at fbr = 6.5^ 
10~'^min~^ (~ 154;t87^ iniii), 7 = 0.3lo2, and /? = 2.11q5 (errors here are 90% confidence 
limits). This solution is plotted over the observed SF in Fig. [2] In Fig. |3] we show the 
confidence contours in the f3 — fbr plane at the best fit value for 7 = 0.3. When the three 
dimensional distribution is marginalized over two parameters, we arrive at the following 
likeliest values and their 90% confidence limits for the individual parameters: fbr 



'•^-5.6 ' 



10 ^min ^ 



77 mmj, 7 



0.6toi, and /3 



9 1 +0.8 



The existence of a break in the power-law PSD of Sgr A* can also be inferred by a 
simple argument: Table 1 reveals that the mean flux of Sgr A* during one night stays 
roughly the same over the 4 years of observations. If the power-law PSD extended unbroken 
over decades, one would expect to observe very different mean fluxes over a time baseline 
of year^ Our detailed MC simulations serve the purpose of quantifying the location of the 
power-law break. The same simple argument also holds true for the mean X-ray flux of 
Sgr A*. 



Our result of a break frequency at ~ 6.5 • 10 ^min" 



^ is strikingly different from the 
prediction of the scaling relation proposed by McHardy et al. (2006). Using a mass of 
M = 4 ■ 10^ Mq (Ghez et al. 2008) and a bolometric luminosity of 2 ■ 10^^ erg/s (Narayan et 



al. 1998), the relation by McHardy et al. (2006) predicts a break frequency at the order of 



10-9 min"^ for Sgr A*. 



Interestingly, if the term with the bolometric luminosity is neglected in the McHardy 



et al. (2006) relation, i.e. the bolometric index (called 'B' in their paper) is set to zero, 
their relation predicts a break timescale of ~ 110 min for Sgr A*. Also, our deduced break 
timescale fits the expected value if linear scaling with BH mass is assumed from the typical 
break timescales observed in the high/soft and low/hard states of Cyg X-1; see Fig. |4j This 
suggests that, while the inclusion of the bolometric luminosity as a free parameter improves 
the fit in a limited luminosity range, it is in fact not the true physical correction factor that 
can explain the scatter around the linear mass scaling. 

It is, however, uncertain how meaningful a comparison of Sgr A*'s break timescale with 
the breaks of the much higher luminosity AGN used to derive the mass-luminosity-timescale 



relation of McHardy et al. (2006) is. These Seyferts are almost certainly in a different 



^Thc high-frequency power-law slope of Sgr A* is close to 2. In this simple case of Brownian motion - 
assuming there is no break in the PSD - the intrinsic standard deviation of the flux scales with the square 
root of time. 
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accretion state. They have a prominent optically thick accretion disk and much higher 



accretion rates. The McHardy et al. (2006) relation seems to work well when extrapolated 



to soft state BHXRBs, which are at comparable Eddington accretion rates to the Seyferts, 
but for the low- accretion rate hard states the comparison is less obvious. In fact, work by 



Gierlinski et al. (2008) suggests that in the hard state the high-frequency shape of the PSD 



is constant despite large changes in luminosity. Therefore, we conclude that the break in the 



NIR PSD of Sgr A* implies that the McHardy et al. (2006) relation with luminosity scaling 
applies to soft state AGN, but that in the hard state there is only mass scaling. Then Sgr A* 
seems to fit in well with the hard state BHXRBs, consistent with its low luminosity. 

Future work will further elucidate the relationship between variability processes across 
the range of BH masses. Our work presented here extends the set of supermassive BHs 
with determined break frequencies. The BH in our Galactic center is especially important, 
as it is the most underluminous source (in terms of Eddington luminosity) accessible to 



observations, and its mass determination is the most precise for supermassive BHs (Ghez et 



al. 2008 Gillessen et al. 2008). 
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Table 1. 


Summary of the data 




Date 


Telescope^ 


Filter Duration 


Mean Flux'' 


published in 


(UT) 






(min) 


(mjy, de-reddened) 




2004 July 06 


Keck 


L' 


165 


2.64 


Hornstein 2007 


2004 July 07 


Keck/VLT 


L'/K 


532 


4.05 


Hornstein 2007/Eckart+ 2006a 


2004 July 08 


VLT 


K 


360 


2.00 


Eckart+ 2006a 


2004 July 26 


Keck 


K 


42 


4.43 


Gliez+ 2008 


2005 July 31 


VLT/Keck 


K 


591 


4.01 


Meyer+ 2008 


2006 May 03 


Keck 


K 


140 


5.53 


Do-I- 2008 


2006 June 20 


Keck 


K 


125 


4.57 


Do-I- 2008 


2006 June 21 


Keck 


K 


164 


3.62 


Do-I- 2008 


2006 July 17 


Keck 


K 


189 


2.86 


Do-I- 2008 


2007 May 18 


Keck 


K 


84 


5.53 


Do-I- 2008 


2007 August 12 


Keck 


K 


57 


3.05 


Do-I- 2008 


2008 May 15 


Keck 


K 


153 


4.57 


unpublished 


2008 June 02 


Keck 


K 


160 


4.00 


unpublished 


2008 July 24 


Keck 


K 


179 


3.51 


unpublished 



Note. — The first three entries are used for the high- to mid-frequency coverage, see Fig. [T] 
''Keck means NIRC2 at Keck II; VLT means NACO at VLT (Yepun). 

''The fluxes are K-band fluxes. When the L' filter was used, we scaled the flux down to the K-band level using 
a spectral index of a = —0.6 | |Hornstein et al.|2007^ . 
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UT hour 



Fig. 1. — The de-reddcncd flux of Sgr A* on 2004 July 06 - 08. The abscissa shows the hours 
elapsed since July 06 00:00 UT. Error bars are omitted for clarity; a sample is shown in the 
upper left corner. It is ±0.75 mJy as determined from the structure function. The two big 
gaps mark the time of daylight. 
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o 




10"^ 10° 10^ 10^ 10*^ 10 

Time Lag [min] 



Fig. 2. — The structure function of Sgr A*. Diamonds represent the high- to mid-frequency 
data shown in Fig. [l| crosses show the mid- to low-frequency data (see Table [l|). The solid 
line is our best fit with Z^. = 6.5 • 10~^min^^, 7 = 0.3, and [3 = 2.1, see Section 4. Please 
note that these parameters describe the PSD and cannot simply be read of the SF depicted 
here. Error bars are determined from our Monte Carlo simulations as described in Section 3. 



10"^ I , , , , I , , , , I , , , , I , , , , I , , , , 

0.5 1.0 1.5 2.0 2.5 3.0 

high — frequency PSD index beta 

Fig. 3. — Confidence contours showing the errors on the best fit parameters. A shoe through 
the three dimensional parameter space at the best fitting 7 = 0.3 is shown. The hnes indicate 
the 99%, 95%, 68%, and 15% rejection probabihty levels. 



I I '-L—l ■ ■ ' I I I ' ' I I I I I I I ' i ' I I I I I I I I I I ' I ' 

^0"^ 0.01 0.1 1 10 100 

PSD break tLme— scale (days) 



Fig. 4. — Sgr A*'s break timescale as reported in this paper over plotted onto Figure 11 of 



Uttley & McHardy (2005), which shows BH mass versus PSD break timescale for various 



AGN. The mass of Sgr A* has been taken from Ghez et al. (2008). Its uncertainty corre- 
sponds to the height of the black square. Filled circles mark masses determined from optical 
reverberation mapping, open circles represent masses determined using other methods. The 
straight lines represent the expected relations if linear mass scaling is assumed from the 
typical timescales observed in the high/soft (dashed line) and low/hard (dotted line) state 
of the BHXRB Cyg X-1 (assuming 10 Mq for its mass). 



